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We present a nonparametric approach to reconstruct the interaction between dark energy and 
dark matter directly from SNIa Union 2.1 data using Gaussian processes, which is a fully Bayesian 
approach for smoothing data. In this method, once the equation of state (w) of dark energy is 
specified, the interaction can be reconstructed as a function of redshift. For the decaying vacuum 
energy case with w = —1, the reconstructed interaction is consistent with the standard ACDM 
model, namely, there is no evidence for the interaction. This also holds for the constant w cases from 
—0.9 to —1.1 and for the Chevallier-Polarski-Linder (CPL) parametrization case. If the equation of 
state deviates obviously from —1, the reconstructed interaction exists at 95% confidence level. This 
shows the degeneracy between the interaction and the equation of state of dark energy when they 


get constraints from the observational data. 


PACS numbers: 95.36.+x, 98.80.-k, 98.80.Es 


I. INTRODUCTION 


It has been more than fifteen years since the universe 
was found in accelerating expansion (al . However, it is 
fair to say that its origin is still not yet clear. A possi- 
ble explanation of this cosmic acceleration is provided by 
the introduction of a fluid with negative pressure called 
dark energy (DE). The simplest dark energy candidate 
is the cosmological constant A with the equation of state 
w = —1. The tiny cosmological constant together with 
the cold dark matter (CDM) (called the ACDM model) 
turned out to be the standard model which fits the cur- 
rent observational data sets consistently. In spite of this 
success, however, it is faced with the fine-tuning prob- 
lem | and the coincidence problem. The former arises 
from the fact that the present-time observed value for 
the vacuum energy density is more than 120 orders of 
magnitude smaller than the naive estimate from quan- 
tum field theory. The later is the question why we live in 
such a special moment that the densities of dark energy 
and dark matter are of the same order. 

Many attempts have been made to tackle those issues, 
including introducing “dynamical” dark energy or mod- 
ifying general relativity at the cosmic scales. In particu- 
lar, to alleviate the coincidence problem, an interaction 
between dark energy and dark matter has been consid- 
ered extensively in the literature, for example, in 
and references therein. However, usually the interaction 
form has to be assumed. The assumed form will lead 
to a bias when observational data are used to give the 
constraint on the interaction. In the authors inves- 
tigated the possible interaction in a way independent of 
specific interacting forms by dividing the whole range of 
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redshift into a few bins and setting the interacting term to 
be a constant in each redshift bin. It was found that the 
interaction is likely to cross the noninteracting line and 
has an oscillation behavior. Recently, Salvatelli et al. 
showed that the null interaction is excluded at 99% con- 
fidence level (C.L.) when they added the redshift-space 
distortions (RSD) data to the Planck data for the decay- 
ing vacuum energy model (a class of interaction of dark 
sectors). They parametrized the interaction term to be 
of the form qHp and also subdivided the redshift into 
four bins with q(z) = qili = 1,...,4). More recently, the 
authors of Ref. have reconstructed the temporal evo- 
lution of the coupling strength between dark matter and 
vacuum energy, a(a) in a nonparametric Bayesian ap- 
proach using the combined observational data sets from 
the cosmic microwave background radiation, supernovae 
and large scale structure. It was found that an evolving 
interaction can remove some of the tensions between dif- 
ferent types of data sets, and is favored at ~ 95% C.L. 
if the baryon acoustic oscillations measurements of the 
BOSS Lyman-a forest sample are included. Thus, it is 
quite interesting to see whether there exists some signa- 
ture of the interaction from the observational data in a 
model-independent way. 

In this paper we present a nonparametric approach 
to reconstruct the interaction term between dark energy 
and dark matter directly from the observational data us- 
ing Gaussian processes (GP). GP is a model independent 
method to smooth the data. We set the nonparameter- 
ized interaction term Q(z) as a function of redshift and 
reconstruct it from SNIa Union 2.1 data sets. We will 
consider three cases, the decaying vacuum energy case 
with w = —1, the wCDM model and the Chevallier- 
Polarski-Linder (CPL) parametrization of dark energy, 
respectively, and pay special attention to the first case as 
discussed in . 

This paper is organized as follows. In Sec. [I] we give 
the interacting models of dark energy and dark matter 
in a flat universe. In Sec. we briefly introduce the 
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Gaussian processes method and create a mock data set 
to demonstrate the reliability of the GP reconstruction 
method. Then we apply it to the real data Union 2.1 in 
the decaying vacuum energy case, followed by different 
values of w for comparison and the CPL case. We give 
some discussions and make conclusions in Sec. 


I. THE MODEL 


In a flat universe with an interaction between dark en- 
ergy and dark matter, the Friedmann equation describing 
the evolution of the universe is given by 


8rG 


H’? = z (Om + ppg), (1) 


where pm denotes the energy density of dark matter and 
ppe the energy density of dark energy. However, the 
conservation equations are changed to be 


Pm or 3H pm = =Q, (2) 


bpe + 3H (1+ w)ppe =Q, (3) 


where w is the equation of state of DE, H is the expan- 
sion rate of the universe, and Q describes the interaction 
between dark matter and dark energy. When Q = 0 and 
w = —1, it recovers the standard ACDM model. Unlike 
most of the ways to parametrize the interaction term Q 
using such a form gH p Bis), here we use Q(z) itself and 
want to reconstruct it directly from data using a model- 
independent method. As assumed in (15), the 4-vector Q” 
is proportional to the 4-velocity of dark matter. Hence 
the perturbation of the interaction is not considered in 
this paper. 

Combining the Friedmann equation and the conserva- 
tion equations (87G = 1 throughout the paper), we can 
obtain 


/ 
-wQ =2 (aH? + H?H" — Znen) +2)? 
W 


3H?w' 
w 


— 26 + 3w)H?H’ — | (1+2) 


+9(1 + w)H?, (4) 
where the prime denotes the derivative with respect to 
redshift z. Note that here we also assume w is time de- 
pendent. For convenience, we use a dimensionless q to 


characterize the interaction, i.e., Q = qH. Note that 
here q is not the usual deceleration. In this case, we have 


f 
—waH? =2 (ai? + H?H" — Znen) (1+2)? 
w 


3H’w' 
w 


— 26 + 3w)H?H' — | (1+ 2) 


+ 9(1 + w)R’. (5) 


The luminosity distances at redshift z can be expressed 
as 


A= of de. (6) 


Writing D(z) = (Ho/c)(1+ z)~1dz(z) as the normalized 
comoving distance, we can arrive at 


3p!” D" wD" 3 
D" 3u" 
+ 26 + 3w) + i (1+2) 


pa (7) 


From this, we see that using the observed distance- 
redshift relationship D(z), one can reconstruct the in- 
teraction, once the equation of state w of dark energy is 
given. 


II. RECONSTRUCTION METHOD 


In order to reconstruct the interaction using current 
data sets, we should find a model-independent method 
to reconstruct D(z) and its derivatives. While there 
are several methods such as principle component anal- 
ysis H2, Gaussian smoothing Æl ba and Gaus- 
sian processes RH22, in this paper we will reconstruct 
D(z) and its derivatives more precisely by using the GP 
method. 


A. Gaussian processes 


The Gaussian processes allows one to reconstruct a 
function from data without assuming a parametrization 
for it. We use Gaussian processes in Python (GaPP) 
to derive our GP results. This GP code has been applied 
in many papers RH33. The distribution over functions 
provided by GP is suitable to describe the observed data. 
At each point z, the reconstructed function f(z) is also 
a Gaussian distribution with a mean value and Gaus- 
sian error. The functions at different points z and Z are 
related by a covariance function k(z, Z), which only de- 
pends on a set of hyperparameters / and a+. Here £ gives 
a measure of the coherence length of the correlation in 
x-direction and of denotes the overall amplitude of the 
correlation in the y-direction. Both of them will be op- 
timized by GP with the observed data set. In contrast 
to actual parameters, GP does not specify the form of 
the reconstructed function. Instead it characterizes the 
typical changes of the function. 

The different choices of the covariance function may 
affect the reconstruction to some extent. The covariance 
function is usually adopted as the squared exponential 
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form 


k(z,Z) =of? ( = a 8 

(z, Z) = of“ exp TE (8) 
But it is not always the suitable choice. Here we take the 
Matérn (v = 9/2) covariance function 


k(z, Z) = of? exp =) 


3lz—2| 27(z-2) 
|z al (z— 2) 


Z 70? 


18|z—2|°? 27(z—2)* 
$ |z — ž| (z a), 


7B | 350 (9) 


according to the analysis made in [29}, where they con- 
sider various assumed models and many realizations of 
mock data sets for a test and conclude that the Matérn 
(v = 9/2) covariance function can lead to more reliable 
results than all others when applying GP to reconstruc- 
tions using D measurements. 

Following Refs. (27, Bo, in which the detailed techni- 
cal description of GP can be found, we reconstruct the 
interaction between dark energy and dark matter using 
the SNIa Union 2.1 data set . Before that we will first 
show the reliability of the GP method. 


x [1+ 


B. Mock data 


To demonstrate the ability of the GP method to distin- 
guish different models and recover the correct behaviors 
of the models, we create mock data sets of future SNIa 
according to the Dark Energy Survey (DES) for two 
fiducial models: the standard ACDM model and a toy de- 
caying vacuum model: ppg = 3aH with w = —1. Here 
we set Qmo = 0.3 for both. 

For the ACDM model, it is straightforward to calcu- 
late the Hubble parameter H/Ho and then obtain the 
simulated data of D(z) easily. Because there is no inter- 
action, we just want to check whether the reconstruction 
with this simulated data can recover q = 0. 

As for the decaying vacuum model, on the one hand we 
should calculate the Hubble parameter and on the other 
hand, we must get q(z) for this fiducial model. The Hub- 
ble parameter is simply H/ Ho = 1—Qmo+2mo(1+z)?/?. 
The fiducial interaction is gHp? = —3(1 — Qmo)(1 + 
z)H'H. The rest of the task is to reconstruct q(z) from 
the Hubble parameter, and then test how well the re- 
constructed q(z) agrees with the fiducial one. If the GP 
method can recover both of the fiducial models and has 
the ability to distinguish them, we can demonstrate that 
GP is a valid method in the reconstruction for our pro- 
pose. In the following we will create the mock data sets. 

The DES is expected to obtain high quality light curves 
for about 4000 SNe Ia from z = 0.05 to z = 1.2 in the 
next five years. From Table 14 in we can calcu- 
late the errors of D: op and the corresponding num- 
bers of SNe Ia for each redshift bin. At every redshift 


point z, D(z)sim is sampled from the normal distribution 
D(Z)sim ~ N(D(Z) pias 7D) where D(Z)sim is the simu- 
lated data of D(z) and D(z) fia is the theoretical value 
from the fiducial model. Thus we create the mock data 
sets for these two fiducial model. For each of them we re- 
construct D(z) and its derivatives and then apply to the 
reconstruction of the interaction term q(z). Note that 
each of the reconstructed D(z) and its derivatives at ev- 
ery redshift point is a distribution with a mean value and 
the error regions. However, D(z) and its derivatives are 
not independent but correlated by a covariance between 
them. Using the reconstructed D(z) and its derivatives 
also the covariance matrix between them, we can apply 
Monte Carlo sampling to determine the q(z) in Eq. 
for a given w at each point z which we want to recon- 
struct. The detailed description of the covariance matrix 
can also be found in Ri. 

We can see from Fig. []that the D(z) and its deriva- 
tives are reconstructed very well from the mock data sets 
assuming the ACDM model. The dashed blue line is the 
mean of the reconstruction and the shaded blue regions 
are the 68% and 95% C.L. of the reconstruction. The er- 
rors for higher derivative of D(z) are a litter larger than 
the lower ones. As expected, at higher redshifts the er- 
rors become large due to the poor quality data in that 
region. 

Since our reconstruction involves D” (z) whose errors 
are even lager in the higher redshift, and will definitely 
lead to somehow uncontrollable large errors in high red- 
shift regions when apply it to the Monte Carlo sampling 
in Eq. (7. For a better show, we introduce a prefac- 
tor (1+ z)~” to q(z), that is, ¢(z) = ¢(z)(14+<2)~”. The 
choice of n is somehow arbitrary, we take n = 6 here. The 
motivation to introduce the prefactor is comprehensible: 
we just focus on the interaction in the low and medium 
redshift range because the quality of observed data in 
the higher redshift range is so poor that it gives a weak 
constraint on the reconstructed interaction. Such a pre- 
factor is just considered as a scale transformation with 
respect to redshift, which does not influence the recon- 
struction of g(z) in the low and medium redshift ranges 
significantly and provides a better show in the higher red- 
shift range. Moreover, our aim is to examine the evidence 
of the interaction, namely, we are testing the equality of 
the quantity with zero. As a consequence, we are free 
to do this without loss of generality like what Ref. [30] 
has done for the null test of the ACDM model. In the 
rest of this paper, we use g as our interaction term to 
test the interaction and we mainly focus on the low and 
medium redshift ranges because of the poor quality of 
the observed data in high redshift ranges which will not 
give a good constraint for the interaction term. 

Figure[2] shows that the reconstructed interaction ¢(z) 
is consistent with the fiducial ACDM model nicely, falling 
in the lo limit. The reconstruction in the redshift range 
from 0 to 0.6 is better than that in the high redshift where 
the error is large. Further, we see from Figs B] and [4] 
that the reconstructions of the decaying vacuum model 
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FIG. 1: Gaussian processes reconstruction of D(z), D'(z) (top), and D” (z), D™ (z) (bottom) obtained from a mock data set of 
future DES and assuming the ACDM model with Qmo = 0.3 (red line). The dashed blue line is the mean of the reconstruction 
and the shaded blue regions are the 68% and 95% C.L. of the reconstruction, respectively. 
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FIG. 2: Reconstruction of ¢(z) (dashed line) from the mock data set of future DES and assuming the ACDM model with 
Omo = 0.3. The shaded blue regions are the 68% and 95% C.L. of the reconstruction. 


also recover the fiducial model very well, and obviously 
deviate from the ACDM model. This shows that the 
Gaussian processes can capture both the two models very 
well and correctly distinguish between them. 


C. Reconstruction using Union 2.1 data 


We now apply the GP method to the real data, Union 
2.1 data sets , which contains 580 SNela data. Let us 


consider 


m—M +5log a —25=5log[(l1+2)D], (10) 
with Ho = 70 km/(s Mpc), following B4. Actually our 
results are not sensitive to the values of Hp. The values 
of D depend on both Ho and the absolute magnitude M. 
We can fix the Hp and only consider the uncertainties in 
M as adopted in Ref. . We transform the distance 
modulus m — M given in the data set to D and set the 
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FIG. 3: Gaussian processes reconstruction of D(z), D'(z) (top), and D” (z), D™ (z) (bottom) obtained from a mock data set of 
future DES and assuming a toy decaying vacuum model: ppg = 3aH with w = —1 and Qmo = 0.3 (green line). The dashed 
blue line is the mean of the reconstructions and the shaded blue regions are the 68% and 95% C.L. of the reconstruction, 


respectively. The ACDM model is also shown (red dotted). 
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FIG. 4: Reconstruction of ¢(z) from a mock data set of future DES and assuming a decaying-vacuum model: ppg = 3aH with 
w = —1 and Qmo = 0.3 (green line). The shaded blue regions are the 68% and 95% C.L. of the reconstruction. The ACDM 


model is also shown (red dotted). 


theoretical initial conditions D(z = 0) = 0 and D'(z = 
0) = 1. First, let us consider the decaying vacuum energy 
case with w = —1 (15). Figure [5]shows the reconstructed 
D(z) and its derivatives from Union 2.1 data. The error 
is a little larger than that for the reconstruction of the 
DES survey due to the smaller number of SNe Ia and 
larger measurement errors. The reconstruction of the 
interaction ¢(z) is shown in Fig.[6] We see that both the 


distance D(z) and interaction ¢(z) are consistent with the 
ACDM model, which implies that there is no evidence for 
the existence of the interaction. 

We now change the equation of state w to see how 
much the differences of w will influence the output of 
our reconstruction. We consider the wCDM model with 
w = —1.006 + 0.045 from the Planck 2015 [36], the 
CPL parametrization w(a) = wo + wa(1 — a) with wo = 
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FIG. 5: Gaussian precesses reconstruction of D(z), D'(z) (top), and D” (z), D" (z) (bottom) obtained from Union 2.1 data sets. 
The shaded blue regions are the 68% and 95% C.L. of the reconstruction. The ACDM model (red line) is also shown. 
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FIG. 6: Reconstruction of g(z) from Union 2.1 data sets. The shaded blue regions are the 68% and 95% C.L. of the reconstruc- 


tion. The red line corresponds to the ACDM model. 


—1.046+9:17 and wa =0.14+5-? from HST Cluster Su- 
pernova Survey 2011 B4 and other constant equation 
of state: w = —0.7, —0.8, —0.9, —1.1, —1.2 and —1.3, 
respectively, for a comparison. All of the results are pre- 
sented in Fig. [7] Here we emphasize that we have con- 
sidered the effects of the errors of w on the reconstructed 
q(z) in the cases of (d) and (e) in Fig. [K] Since the con- 
stant w for the wCDM model from Planck 2015 is very 
close to —1, so the reconstruction is almost the same as 
the w = —1 case. While for the CPL case, the error of 
the reconstruction is a little larger than those in other 


cases because of the poor constraints on wo and wa. We 
see that if w lies between —0.9 and —1.1, G(z) = 0 is 
captured within the 95% confidence region of the recon- 
struction, while the interaction is shown up if the equa- 
tion of state for dark energy deviates significantly from 
—1. This shows the fact that there is a degeneracy be- 
tween the interaction and the equation of state of dark 
energy, as indicated in Eq. (7). 
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FIG. 7: Reconstructions of q(z) for the wCDM model and CPL parametrization. (a) w = —0.7; (b) w = —0.8; (c) w 
(d) w = —1.006 +0.045 from Planck 2015; (e) w(a) = wot+wa(1—a) with wo = —1.046+9:119 and wa = 0.14+9:80; (f) w 
(g) w = —1.2; (h) w = —1.3; The shaded blue regions are the 68% and 95% C.L. for the reconstruction. 
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IV. DISCUSSIONS AND CONCLUSIONS 


We have presented an approach to reconstruct the in- 
teraction between dark energy and dark matter by using 
Gaussian processes. To check the reliability of the GP 
method, we create mock data for two fiducial models, 
one is the ACDM model which has no interaction be- 
tween dark matter and dark energy, the other is a toy 
decaying vacuum energy model with ppg = 3aH, which 
is of the interaction. It shows that the reconstruction 
method by using the Gaussian process works well and 
can capture the features of these two models. 

We then applied the method to the real data from 
Union 2.1 data sets and reconstructed the distance and 
its derivatives, and then the interaction. It was found 
that for the decaying vacuum energy model with w = 
=l E there is no evidence for the existence of the in- 
teraction, namely, the ACDM model is consistent with 
the Union 2.1 data sets within 1ø limits. 

To check the influence of the equation of state of dark 
energy on the method, we have also considered several 
constant values of w from —0.7 to —1.3. The results 
show that g = 0 falls in 95% C.L. of the reconstruction if 
w lies between —0.9 and —1.1. The wCDM model with 
w = —1.006 + 0.045 from Planck 2015 falls in this range 
and q = 0 is within the 1o limits. This also holds for the 
CPL parametrization case. However, as we can see from 
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Figure [7| if w deviates obviously from —1, the interac- 
tion exists beyond 20 C.L.. This reflects the degeneracy 
between the interaction and the equation of state of dark 
energy. 

Note that in our reconstruction method, only the ob- 
servational data on the expansion history of the universe 
can be used. For example, some measurements of Hubble 
parameter can also be combined into the above recon- 
struction. It is certainly of interest to develop a recon- 
struction method for the interaction between dark mat- 
ter and dark energy, in which some kinetic data of the 
universe, for example, the growth factor of large scale 
structure, can be used for this aim. Our paper shows 
that the existence of the interaction between dark mat- 
ter and dark energy found in is mainly due to the 
data of redshift space distortion. 
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